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Unusual electronic property of a Weyl semi-metallic nanowire is revealed. Its band dispersion 
exhibits multiple subbands of partially flat dispersion, originating from the Fermi arc states. Re- 
markably, the lowest energy flat subbands bear a finite size energy gap, implying that electrons in 
the Fermi arc surface states are susceptible of the spin Berry phase. This is shown to be a conse- 
quence of spin-to-surface locking in the surface electronic states. We verify this behavior and the 
existence of spin Berry phase in the low-energy effective theory of Fermi arc surface states on a 
cylindrical nanowire by deriving the latter from a bulk Weyl Hamiltonian. We point out that in any 
surface state exhibiting a spin Berry phase n, a zero-energy bound state is formed along a magnetic 
flux tube of strength $o/2 = he/ (2e). This effect is highlighted in a surfaceless bulk system pierced 
by a dislocation line, which shows a ID chiral mode along the dislocation line. 



I. INTRODUCTION 

It has been proposed recently that a three-dimensional 
(3D) Weyl semimetal phase is likely realized in pyrochlore 
iridates, A 2 Ii"207 (A=Y, or a lanthanide element) P The 
Weyl semimetal has a unique band structure character- 
ized by a set of discrete linearly-dispersive point nodes, 
the Weyl pointsPEl Such a band structure of Weyl 
semimetal is naturally reminiscent of that of graphene, 
its 2D counter part, but unlike 2D Dirac cones: H^d = 
VxOx +P y o-y, 3D Weyl points: H 3D = p x a x +p y <T y +p z o~ z 
are more stable objects; they cannot be trivially gapped 
out by a perturbation (cf. Hid is gapped by a local mass 
term, H' = ma z ). 

Unlike topological insulators also, the Weyl semimetal 
is gapless in the bulk, but when in contact with the vac- 
uum it exhibits a peculiar surface state, somewhat remi- 
niscent of a more established helical surface state of the 
topological insulator. The two surface states are indeed 
both of topological origin, and in this sense this analogy 
between the Weyl semimetal and the topological insu- 
lator is not superficial at all. In the bulk (without a 
surface) the Weyl semimetal is already gapless, but the 
conduction and the valence bands touch only at discrete 
points (Weyl points) in the Brillouin zone - let us assume 
hereafter that there exists only a pair of such point nodes 
for simplicity. In the presence of a surface, an additional 
state appears, localized on the surface, and "enveloping" 
the two point nodes. If one considers the E = cross sec- 
tion of the energy spectrum, this additional surface state 
appears as a line, not necessarily straight but al ways con- 
necting the two Weyl points, forms a Fermi arc! 1 * 2 * 5 * ^ As 
is typically the case with the helical Dirac cone surface 
state of a topological insulator, the existence of this Fermi 
arc envelope state is topologically protected by a bulk 
topological invariant through the so-called bulk/surface 
correspondence. 

As we mentioned earlier, the Weyl semimetal can be 
regarded as a 3D version of graphene. Such an anal- 
ogy in the low-energy electronic property of the hulk is 
naturally extended to that of the surface. The Fermi 



arc state is indeed shown to be a precise 3D analogue of 
the edge states of a graphene nano-ribbon in the zigzag 
edge geometryP These two examples constitute prototyp- 
ical classes of the topologically non-trivial gapless states, 
which are counterparts of the gapped topological insula- 
tor and superconductors, the latter known to be classified 
into the form of a periodic table (ten- fold way) in terms of 
their symmetry and dimension ! 10 ! 11 ! The idea of charac- 
terizing topologically non-trivial gapless or nodal states 
in terms of the topological invariants has bee n int roduced 
and extensively used in the study of 3 He-A.l2E2Hm More 
recently it has been applied to the description of topolog- 
ically non-trivial nodal superconductors G2HZD A periodic 
table analogous to the one used for classifying various 
classes of topological insulators and superconductors has 
been also proposed for those classes of topologically non- 
trivial gapless statesP^l 

The edge/surface state of a topological insulator is of- 
ten referred to be "helical" , indicating that its spin direc- 
tion is locked with respect to its propagating direction. 
Here, in this work we focus on still another unique prop- 
erty of such a helical surface state, i.e., the existence of 
the spin Berry phase. The electronic spin in the sur- 
face helical state shows also spin-to-surface lockingP^HUl 
i.e., the spin is locked in-plane to the tangential surface of 
the real space geometry (e.g., on a cylindrical surface). A 
mathematical description of this spin-to-surface locking, 
the spin Berry phase primarily manifests in the finite-size 
energy gap associated with the surface helical statesP^In 
an infinitely large (or doubly periodic) slab geometry, the 
finite-size energy gap of the surface state due to a finite 
thickness of the slab decays exponentially as a function 
of the thickness. In the case of a rectangular/cylindrical 
nanowire, i.e., when the width of the slab becomes fi- 
nite, and the slab acquires side surfaces, this is no longer 
the case. The phase information of the electronic wave 
function on one surface can be transmitted to that of 
the opposite surface via the (gapless) side surface states. 
Then, the finite-size energy gap decays only algebraically 
as a function of the thickness. The spin Berry phase 
replaces the periodic boundary condition applied to the 
electronic motion around the cylinder by an anti-periodic 
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FIG. 1. Multiple subbands, originating from a pair of Fermi 
arc states connecting the two Weyl points (k z — ±fco = ±7r/2) 
via the Brillouin zone boundary. The subbands have a flat 
band dispersion, filling the bulk energy gap at k z > ko. A = 
B — t z — 1. The simulation is done for a system of square 
lattice of size: (N x ,N y ) — (24,24), whereas periodic in the 
z-direction. The inset shows details of the spectrum around 
k z — ko and E — 0. 



boundary condition, leading to half-integer quantization 
of the orbital angular momentum around the cylinder.^ 
Such an energy-gap due to phase coherent motion of an 
electron around the (cylindrical) surface is sensitive to 
introduction of a 7r-flux tube piercing the nanowire. The 
sensitivity to 7r-flux is a fingerprint of the existence of 
spin Berry phase, which might be directly triggered ex- 
perimentally in an Aharonov-Bohm type measurement 
recently performed in a system of topological insulator 
nanowire. 1 ^ 

As mentioned earlier, the Fermi arc states resemble the 
edge states of a zigzag graphene nano-ribbon from the 
viewpoint of topological classification. But still, as we 
demonstrate in this paper, the Fermi arc states exhibit, 
unlike a 2D graphene layer, specific spin Berry phase. In 
this regard the Fermi arc states show a stronger resem- 
blance to the helical surface states of a 3D topological in- 
sulator. This paper reveals smoking-gun features of such 
spin Berry phase in the Fermi arc states. We first demon- 
strate in Sec. II that in a nanowire geometry electrons in 
the Fermi arc surface states show multiple subbands of a 
partially flat dispersion (see FIG. I), but they are suscep- 
tible to a finite-size energy gap associated with the spin 
Berry phase. We then confirm in Sec. Ill the existence 
of this spin Berry phase in the surface effective Hamilto- 
nian by deriving it from the bulk effective Hamiltonian. 
In Sees. IV and V we analyze the system's response to 
an Aharonov-Bohm flux as well as to introduction of a 
screw dislocation, confirming the existence of spin Berry 
phase. We set H = 1, unless otherwise mentioned. 



II. MULTIPLE FLAT SUBBANDS 

Let us consider a simple model of a Weyl semimetal 
with a single pair of Weyl nodes on the A^-axis!^ 

H = A(k x <T x + k y <Ty) + M(k)a z , (1) 

where k — (k x , k y , k z ), 

M(k) = A(k z ) + B(k 2 x + k 2 y ), (2) 

and we choose 

A(k z ) = 2£ z (cos k z — cos ko). (3) 

This is a long- wavelength effective Hamiltonian regarding 
the motion in the x- and y-directions, whereas it can be 
regarded, as for the motion in the z-direction, a tight- 
binding Hamiltonian. Or, by making the replacements: 

k x y sin k x , 

ky v sin ky , 

k 2 x -> 2(1 - cosfer), 

k 2 y -> 2(1 - cosfc y ), (4) 

the same model can be viewed as a 3D square-lattice 
tight-binding model. We will employ this square-lattice 
implementation for numerical simulations. 

The energy spectrum of this effective Weyl model is 
characterized by a pair of Weyl points at k = (0, 0, ifcrj), 
exhibiting a conic dispersion a round them. Besides, a 
pair of Fermi arc states appear GEMS when we put this 
system into a slab, say, bounded by the two surfaces, one 
at x — and the other as i = L x , parallel to the z- 
axis. Let us fix the parameters such that t z > 0, B > 0. 
Then for -k < k z < k , A(k z ) > 0, i.e., A(k z )/B > 0. 
This means that a cross section of Eq. |T]) at a fixed k z 
in the above range, describes a trivial 2D band insula- 
tor. Whereas, for ko < k z < n and — ir < k z < —ko, 
A(k z ) < 0, i.e., A(k z )/B < 0. Then, a similar cross sec- 
tion of Eq. ([I]) at k z in one of these ranges, describes a 
topological (quantized anomalous Hall) insulator with a 
chiral edge mode. There appears one single chiral branch 
on the x = side, and another on the x = L x side, 
propagating in the opposite directions; — y and +y, re- 
spectively. These chiral modes show a linear dispersion, 
therefore, of opposite sign, and cross at k y — and at 
E = 0, forming a "X-shaped" dispersion E = E(k y ). If 
one allows k z to vary continuously, then the locus of such 
an X-shapcd dispersion forms two planar membranes in 
the (k y , k z , i?)-space, always crossing at k y = (on the 
fc z -axis) and at E = 0. The locus of the crossing point is 
the Fermi arc, connecting the two Weyl points k z = ko 
and k z — — ko via the zone boundary. Both ends of the 
two planar membranes are closed by half-conic structures 
which appear "beyond" the Weyl points; k z < ko and 
— ko < k z . The entire manifold thus formed envelops the 
two Weyl cone regions. 

Let us then further restrict the system into a nanowire 
geometry; the system is restricted not only between x = 
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and x = L x but also between y = and y = L y . The 
energy spectrum E = E(k z ) of such a Weyl semi-metallic 
nanowire is shown in FIG. [I] The spectrum shows a series 
of flat subbands: 



E{k z 



(5) 



which are remnant of the two planar regions of the Fermi 
arc manifold. The flatness of the subbands stems from 
the fact that the membrane state has no dispersion in 
the fc z -direction. These multiple subbands form circular 
chiral modes, carrying a spontaneous persistent current 
around the surface of the wire. Note that such circu- 
lar chiral modes do not appear if one considers a wire 
perpendicular to the z-axis. In that case Fermi arc type 
surface states do appear on the surface parallel to the 
z-axis, but they disappear on the side normal to the z- 
axis. Therefore, the surface states cannot completely 
wrap the wire. This is quite contrasting to the topo- 
logical insulator surface states. The latter, protected by 
the "strong" bulk/edge correspondence, appear, irrespec- 
tively of the shape and direction of the surface, and con- 
sequently cover the entire surface. Here, the Fermi arc 
subbands are indeed in one-to-one correspondence with 
the structure of Weyl points in the bulk spectrum. This 
correspondence is, however, "weak" in the sense that it 
depends on the direction of the surface. As we mentioned 
earlier, both the existence and flatness of the Fermi arc 
subbands are also topologically protected. 

What might be counter-intuitive in FIG. [T] is that 
the lowest conduction and the highest valence subbands, 
E±i(k z ) are still separated by a finite size energy gap (see 
its inset). We see later (FIG. [5] upper panel) that one 
can actually close this gap by introducing a flux w pene- 
trating the cylinder. In the light of our knowledge on the 
helical surface states of a 3D topological insulatorpSUll 
such a behavior may be naturally attributed to the exis- 
tence of spin Berry phase. 

Let us come back to FIG. [T] and look into the sub- 
band spectra. On the k z < ko side, the behavior of 
E(k z ) directly results from the bulk spectrum. Indeed, 
for A{k z )/B > 0, only bulk solutions are possible. The 
"bulk" solutions are the solutions of Eqs. ([I]), (|2| and 
([3]) satisfying the boundary condition: 



■0(0, y) = ip(L x , y) = ip(x, 0) = ip(x, L y ) 



(6) 



The "plane- wave" solutions of Eqs. Q, ^ and ^ are 

|fc±) planc = e^ x+k ^\k z )\d{k)±), (7) 

where fc = (k x ,k y ,k z ), and \k z ) is a Bloch state with a 
crystal momentum k z in the z-direction. |d(fc)±) repre- 
sents a spin state pointing in the direction of d{k) = 
(Ak x , Ak y , M (fc)) , either parallel or anti-parallel, de- 
pending on the index ±. The plane- wave solution Q 
has an energy eigenvalue, 




E(k) = ±JA 2 (k x + kl) + M (fc) 2 



(8) 



FIG. 2. Behavior of multiple flat subbands in the presence 
of a flux tube <E> = <E>o/2 = hc/(2e) [$o = hc/e] penetrating 
the nanowire. Upper: the total flux tube $ = <3>o/2 pierces 
a single plaquette at each layer. Lower: the flux tube is split 
into two, each carrying $/2 = $o/4 in order to avoid a zero- 
energy state bound to the flux tube. N x = N y = 24. 



For large enough A(k z ) such that only the third com- 
ponent of d(k) dominates, i.e., d(k) \\ z, and one can 
satisfy the boundary condition Eq. (|6| by a simple su- 
perposition of |fc±) with fc = (±fe x , iky, k z \ i.e., 

|fc±)buik ^ sm(k x x) sin(fcj,y)|fe^)|z±), (9) 

where fc = (n x n/L x , n y Tt/L y , k z ) with n x , n y = 1, 2, ■ ■ • . 
Lowest bulk subbands correspond to (n x ,n y ) = (1,1), 
(n x ,n y ) = (1,2) and (2,1), {n x ,n y ) = (2,2), (n x ,n y ) = 
(1,3) and (3,1), etc. In the crossover regime k z ~ k , 
this simple picture is no longer valid, but the subbands 
may be still classified by these quantum numbers. 

On the k z > ko side, the Fermi arc subbands ([5| appear 
in addition to these bulk solutions. The lowest energy 
subband solutions in the bulk merge into the Fermi arc 
(surface) subbands in the crossover regime. The Fermi 
arc subband solutions are solutions of Eqs. 0, ([2| and 

twith ([6| with k x and k y being a complex number, 
eir wave functions are localized on the surface of the 
wire. Last but not the least, the Fermi arc subbands 
([5]) show a flat spectrum, which appear below the bulk 
solutions; \E\ < \A(k z )\, and exist only in the regime: 
A(k z )/B < 0. 
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In the next section, we investigate the nature of such 
surface Fermi arc solutions. We derive a low-energy effec- 
tive Hamiltonian which involves only the surface states. 
It will become clear that the Fermi arc solutions indeed 
emerge from the bulk effective Hamiltonian, Eqs. (JlJ, 
([2]) and p|, but importantly, with the spin Berry phase, 
which explains the finite size energy gap which we have 
seen in FIG.[U 



III. DERIVATION OF THE SPIN BERRY 
PHASE — THE SURFACE EFFECTIVE 
HAMILTONIAN 

To clarify the nature of spin Berry phase, we consider 
here a cylindrical nanowire of radius R extended along 



the z-axis: 



y < R. We start from the same bulk 



effective Hamiltonian, Eqs. 0, ([2| and ([3]), but in or- 
der to extract relevant information on the surfa ce ele c- 
tronic states we divide it into two componentsp^^HUl 
H = H± + H\\, where H± (H\\) describes electronic mo- 
tion perpendicular (tangential) to the cylindrical surface. 
Eqs. ([!]), |2| and ^ represent an effective theory for 
k x ,k y <C 1, but there is no restriction on k z . Here, we 
consider the case: k z > ko such that A(k z )/B < 0, and 



expand it as k z 



p z . Introducing, 



k<p — 



dr ' 



krh = 



.1 d 



(10) 



conjugate to the cylindrical coordinates: r 
4> = arctan ^, one can express H± and H\i as 

H ± = Hj_(k r ) = H\, , ,(o), 



H {{ =H ]l (k^, Pz ). 



(11) 



In order to derive the surface effective Hamiltonian, we 
first have to construct a base solution, the Fermi arc so- 
lution in the present case, satisfying the given boundary 
condition, 



\iP(r = R,<f>,z)) = 



(12) 



Such a base solution is found by solving the electron dy- 
namics perpendicular to the surface, 



where H± reads explicitly 



H 



M ± Ak r e~ 1 ' 
Ak r e^ -M x 



Here, we have decomposed the mass term into 

M x ~ A(k^) + Bk 2 r , 

M11 = Bki - 2t z sin kf ] p z . 



(13) 



(14) 



(15) 



The Laplacian in the cylindrical coordinates has another 
contribution, (l/r)d/dr. Here, we neglect this first-order 



derivative term, keeping only the term <9 2 /<9r 2 .^ This is 
justified when the radius R of the cylinder is suf ficien tly 
larger than the penetration depth see Eq. (16)] of 

the surface state. 

We search for solutions of Eq. (JT3J) , which has an en- 
ergy Ex, in the range: —A(k z ) < E± < A(k z ), and takes 
the following form, 



IV>j 



e K(r-R)\ E 



where 



\E±,k) = 



±,K) 



(16) 



(17) 



and k > 0. For a given energy E±, k has two positive 
solutions, k = k±, satisfying 



El 
M x (k) 



A-Bk 2 



2„2 



A K 



(18) 



Composing a linear combination of these two base solu- 
tions, one can construct a hypothetical wave function, 

IVm) = Cl e^ r -V\E ± ,K+) +c 2 e*-( r -V\E ± ,K_) 7 (19) 

which should be matched with the boundary condition 



(12), i.e. 

det 



E ± +M ± (k + ) E± + M±(k-) 
-iK+Ae 1 ^ -in-Ae** 



= 0. 



(20) 



Note that the two wave functions \ip) and \i>i_) are related 
by \tp) — V'lK'/') z)\if}±)- Since k + 7^ the condition (20) 
simplifies (after some algebra) to 

E ± {A +E X ) = 0. (21) 

Recall that at fixed k z the Fermi arc solution appears 



in the bulk gap: -A(k z ) < Ej_ < A(k z ), i.e., Eq. ( |21| ) 
imposes E± — 0. Substituting this back to Eq. (18), one 



finds Mj_ = ±kA. For Ej_ = and Mj_ = kA, Eq. ((17|) 
becomes 



\E t 



0, k) = M± 



1 

-ie l 



and the two solutions for n becomes 

A ± y/A 2 + 4BA 



2B 



(22) 



(23) 



which is concistent with the condition: k± > (recall 
that BA < 0). The other choice, M± = —kA is not 
compatible with this requirement. 

Thus the normalized Fermi arc base solution is found 
to be 



\ip x ) = p{r) (e K +( r - fl > - e*"^-^) 



1 



(24) 
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where 
p{r) 



+ k-) e K + {r - R ^ - e K -< r -- R ) 



2nR 



(25) 



(k±R 3> 1 assumed). Eq. (24) is a remarkable result, 
indicating that the surface spin state is 



10-) = 



1 

71 



e -i<t>/2 



(26) 



diagonalizing a spin operator in the direction of with 
an eigenvalue —ft/2, where <f> is a unit vector pointing to 
the azimuthal direction, 



4> = 



— sin < 

cos <b 



(27) 



The electronic spin in the Fermi arc state is locked in 
the direction (anti-) parallel to that of <f>, and when an 
electron goes around the cylinder in the anti-clockwise 
direction, it also rotates, following the curved surface of 
the cylinder, locked in the direction of —(f). After a com- 
plete 2-7T rotation, the electron goes back to the original 
position on the cylinder, the spin also comes back its 
original state, but with an additional phase of tt. This 
may not be clear from Eq. |24|, since it is written in 
the single- valued representation. 1 ^ Yet, information on 
the double-valuedness of spin is safely encoded in the 
surface effective Hamiltonian in the form of spin Berry 
phase, as we see below. Notice also that here, in con- 
trast to t he ca se of helical surface states of a topological 
insulator) 35 ! 36 ! the surface spin state is not helical. It is 
rather "chiral" , pointing to the azimuthal direction of the 
cylinder independently of the value of k z . The spin direc- 
tion is locked indeed anti-parallel to the group velocity of 
the surface mode. Let us finally see such chiral spin-to- 
surface locking leads, indeed, to the appearance of spin 
Berry phase. Using the base solution (24), we calculate 
the expectation value of 



#n = 



-iAe^k, 



-iAe-Wkt 

-Mn 



to find 



tfsurf = (1/>X\HM ± ) = 



(28) 



(29) 



The low-energy electron dynamics along the surface is 
thus determined by the eigenvalue equation, 



H sur { V'H (0, z ) = E \\ i>\\ {4>, z) 



(30) 



where if>\\ (<f), z) = e lk< t , ^e' lk ' z . The periodic boundary con- 
dition around the wire, 



M</> + 2n,z)=M<l),z), 



(31) 



requires that fc^ be an integer. The spin Berry phase 
term, i.e., factor 1/2 in Eq. (29 1 plays, then the role of 



shifting the surface electron spectrum by a half of the 
finite-size energy gap, 



A 
R 



Av 



(32) 



Recall that the origin of this 1/2 spin Berry phase term 
is that the spin in the Fermi arc surface state is locked in 



the direction of Eq. ( 26 1. This is quite contrasting to the 
spin state of the bulk solution; see Eqs. Q and (J9|. In 
the crossover regime k z ~ fco, the bulk spin state \d±) on 
the k z < fco side evolves into the locked surface spin state, 
Eq. (26 1. In parallel with this evolution in spin space, 



the wave function of lowest energy bulk subbands tend 
to become localized around the boundary. In spectrum, 
they merge into the Fermi arc (surface) subbands. 



IV. NUMERICAL CONFIRMATION OF THE 
SPIN BERRY PHASE 

Let us verify the existence of spin Berry phase in nu- 
merical experiments. We reconsider the nanowire geom- 
etry with a square cross section, which has already ap- 
peared in Sec. II. Here, to confirm the existence of spin 
Berry phase, we introduce a flux tube $ piercing the 
nanowire. In FIG.[2]the energy spectrum in the presence 
of an infinitely thin 7r-flux tube, carrying a magnetic flux 
$ = $o/2 = hc/(2e), with $ = hc/e being the unit flux 
quantum, is shown, to be compared with the spectrum in 
the absence of flux (FIG. [TJ. The upper panel of FIG. [| 
shows the spectrum when the total flux $0/2 = hc/(2e) 
pierces, at each cross section, a single plaquette. The 
obtained spectrum shows a gapless dispersion with dou- 
bly degenerate zero-energy states. This behavior is in- 
deed quite contrasting to the gapped spectrum of FIG. 
[l] which is reminiscent of an analogous behavior due to 
the effect of spin Berry phase on the surface of a topo- 
logical insulator.-?' Still it is not simply explained by the 
effective surface theory, Eqs. (29 1, (301, (31) and (32), 



which we have derived in the last section. Notice that in 
the presence of a flux $o/2 penetrating the cylinder, the 
surface effective Hamiltonian ( 29 1 is replaced by, 



A 
R 



d 



resulting in a shift of the spectrum, 



iw($) 



A 
R 



k<h + 



$ 



$0 



(33) 



(34) 



Thus, according to the effective surface electron dynam- 
ics, Eqs. (33), (30), (31) and (34), the subband spectra 



are uniformly shifted both for fc^ positive and negative 
(or null) integer. 

Then, where does the degenerate E = pair come 
from? The answer to this question is almost obvious if 
one looks into spatial distribution of the corresponding 



() 



a 



a 




kz / n 



FIG. 3. Multiple subbands deformed by a crystal dislocation 
with a Burgers vector b = (0,0,6), where 6=1 for the top, 
6 = 2 for the central, and 6 = 4 for the bottom panel, fco = 



tt/3, N x = N y 



16. 



wave function. As a general consequence of the 1 /2 spin 
Berry phase term, a series of plaquettes (aligned in the z 
direction) each penetrated by a 7r-flux tube always hosts 
a zero-energy bound state. (A possibly related effect in 
a strong topological insulator can be found in RefP.) 
The existence of such a bound state might be clear from 
Eq. ( 34 ) . As a crude approximation, one can regard the 



series of plaquettes which accommodates the 7r-flux as a 
cylinder of a radius tq ~ 1 (the lattice constant) . Taking 
into account (though this is irrelevant to the discussion 
here) that the surface state localized around this cylinder 
will have an opposite chirality (propagating direction), 



the effective surface Hamiltonian ( 29 ) , may be modified 
to describe such a bound state as 



bound 



A 



•JL 



(35) 



In Eqs. (33 1, (35) the electrons feel the same flux, only 
the propagating direction is opposite. We may write 
the corresponding wave function as ip\\(4>, z) = e " ll *e lfczZ , 
with a quantum number n associated with the orbital 
motion around the flux tube, rather than to make a 



distinction between the two. Then, Eq. (34) becomes 



E, 



bound 



A>) 



A 



4> 



(36) 



In any case, cancellation of the 1/2 spin Berry phase term 
by the 7r-flux implies the existence of a zero-energy bound 
state. Of course, since ro ~ 1 <C R, in the spectrum of 
Eq. (36) only the E — (n — 0) state is relevant in 
the energy scale A/R of the finite-size energy gap; cf. 
Eq. (34), and appears in the window of bulk energy gap. 



Such a bound state along the flux tube is degenerate with 
the E = subband state of Eq. ( 32 ) with = and 
explain the two-fold degeneracy of E = state in FIG. 
® 

However, if one's purpose is to see simply the effects 
of spin Berry phase, one can avoid this complexity. The 
lower panel of FIG. [2] shows a spectrum when the system 
is always pierced by a tt-Aux tube, with a magnetic flux 

= $0/2 = hc/(2e), but split into two; each of the half 
flux <f>/2 = $o/4 = hc/(Ae) pierces a different plaquette. 
One can still assume a bound state along such a half flux 
tube, and estimate its energy. For a cylinder penetrated 
by a half flux tube, Eq. (36) modifies to, 



E 



bound 



($) 



E 



A 



bound 



($/2) = 



1 $/2 

H f— 

2 $0 

(37) 

Clearly, for $ = $o/2 and ro ~ 1 -C R, there exists no 
bound state in the scale of finite-size energy gap A/R. 
The n = bound state is sent to the high-energy spec- 
trum (one can actually see this in the lower panel of 
FIG. [2]). Thus, the low-energy (at the energy scale of 
A/R) spectrum of such a system is simply determined 
by Eq. ([34]). The obtained data depicted in the lower 
panel of FIG. [2] indeed show such a behavior consistent 
with Eq. pi] ), with a single subband state located pre- 
cisely at E = (see the inset). 

An alternative way to verify the existence of spin Berry 
phase is to see system's response to crystal deformation; 
introduction of a screw dislocation. A screw dislocation 
plays, fundamentally, a role similar to the magnetic flux 
we have considered above but its effect on the sub- 
band spectrum is superficially much different, as shown in 
FIG.[3j Suppose that the underlying crystal is deformed 
by a screw dislocation along the axis of the wire (z-axis); 
its Burgers vector is b = (0, 0, b) (b = ±1, ±2, • • • ). As 
opposed to a magnetic flux which twists the phase of an 



7 



electronic wave function uniformly, a crystal dislocation 
introduces a phase shift which depends on the crystal 
momentum k z , in the direction of the Burgers vector, 



V^ll (<f>, z) 



(38) 



introducing a finite slope into the subband spectrum, 



, , A ( 1 k z b 
E sul ,(k z b) = -lk,+ 2 - — 



(39) 



Notice that the fc^-th subband intersects with the E = 
line at 



k z b = (2fe^ + 1)tt, 



(40) 



similarly to the 7r-flux case, but here this occurs only at 
such discrete values of k z . 

Is that all that a dislocation line does to the subband 
spectrum? No, of course not. In each panel of FIG.[3]one 
can recognize an isolated mode which has a slope opposite 
to all the other subband states. This is again due to a 
bound state formed along a dislocation line. Similarly to 
the case of a tt flux tube piercing a single plaquette at 
each (x, y)-layer, a series of plaquettes penetrated by the 
dislocation line may be regarded as a cylinder of radius 
ro ~ 1. The spectrum of subband states associated with 
such a dislocation line reads 



Ebound(k z b) — 



A 

ro 



kzb 
2ir 



(41) 



Again, for ro ~ 1 only at most a few subbands, satisfying 
a zero-energy condition, 



k z b = (2n + l)7r, 



(42) 



are visible in the relatively small window of the bulk spec- 
trum. Such subbands has a steep positive slope as a func- 
tion of k z opposite to all the other subbands described 
by Eq. (39 1 since R >■ 1; i.e., the bound state(s) along 
the dislocation line is (are) propagating modes. Notice 
that in the two upper panels of FIG. [3] (cases of b = 1,2) 
crossing of the two subbands, ( [39| and ( |4"T] ), occurs at 
k$ = n — 0, whereas in the bottom panel (6 = 4 case), 
the same crossing which falls on the Fermi arc region is 
at k<j, = n = 1, indicating that k z — (3/4)7r. 

We have seen that a (square cross sectional) cylinder 
pierced by a dislocation line should be regarded as a pair 
of concentric cylinders with radii, tq ~ 1 and R 1. 
This allows us to interpret the spectra shown in FIG. [3] 
as a superposition of two contributions; one from the sub- 
bands localized in the vicinity of the outer surface: Eq. 
(39), and the other from a bound state along a disloca- 
tion line: Eq. (41 ). To double check the validity of such 



an interpretation in terms of the bound state along a dis- 
location line, we consider in the next section an extreme 
example, in which only the bound states appear. 




FIG. 4. ID chiral modes along a pair of dislocation lines. 
6=1 (top), b — 2 (central) and 6 = 4 (bottom), ko = 7r/3, 



N x = N y = 24. 



V. PROTECTED ID CHIRAL MODE ALONG A 
DISLOCATION LINE 



Let us finally consider a slightly different geometry in 
which a pair of screw dislocations, one parallel, the other 
anti-parallel with the z-axis, penetrate a triply periodic 
(surfaceless) system; the two dislocation lines are spa- 
tially well-separated. On a planar region bounded by the 
two dislocation lines the crystal is dislocated in the z- 
direction by b. This attributes to each dislocation line a 
pair of Burgers vectors, b — (0,0,6) and b = (0,0,-6). 
As we have already seen, such a situation is readily de- 
scribed by the cylinder model we have considered in Sec. 



IV. A minor but not unimportant difference from the pre- 
vious case is that here the two cylinders are parallel, and 
not concentric. Therefore, on the surface of the two cylin- 
ders the orbital motion of an electron around the cylinder 
is in the same anti-clockwise direction. The low-energy 
electron dynamics on the surface of the two cylinders are 



described by the same surface effective Hamiltonian ( 29 1 . 
Only the Burgers vector differs, and modulates the phase 
of the electronic wave function in different ways; along 
each of the dislocation line, Eq. ( 38 ) should be replaced 
by, 



i>l((f>,z) = e i(«i-fc*&/(27r))0 e ife^ 
</> 2 (<M = e *(«2+fc,fc/(2?r))0 e ^ 

The corresponding bound state spectra read, 



E x {k z b) 



A 
A 



ni 



E 2 (k z b) — I T12 , , 



kzb 

k z b 
2n 



(43) 
(44) 



(45) 
(46) 



Again, since r ~ 1 only the lowest energy subbands, 
satisfying the zero-energy condition, 



k z b = (2m + l)vr, -k z b = {2n 2 + 1)tt, 



(47) 



are relevant in the spectrum. Some concrete examples of 
calculated spectra for such a system are shown in FIG. 
[i] for b = 1 (top), 6 = 2 (central) and 6 = 4 (bottom). 
Each spectrum exhibits a pair of chiral modes, which are 



identified as the states represented by Eqs. ( 43 ) and ( 44 ) 



with n\ and n 2 satisfying Eqs. (47) in the Fermi arc re- 
gion: ko < k z < 7r. For the two upper panels (cases of 
6 = 1,2), 774 = 0, ri2 = — 1) i-c, the two chiral modes 
intersect at k z — ir/b and at E — 0. As for the last 
panel (6 = 4), ri\ = 1, n 2 — —2, indicating that the in- 
tersection occurs at k z = (3/4)7r. Note that this type of 
gap closing always occurs at E = and at the same k z 
points uniquely determined by the Burgers vector. Such 
a feature is model independent, and in this sense these 
chiral modes are protected. Notice, in contrast, that pro- 
jection of the 3D Weyl point onto the ID Brillouin zone 
(k z = 7r/3 in FIG. [4]) is gapped by the screw dislocation. 

It is also interesting that such a pair of zero energy 
bound states have a dispersion in the fc z -direction; they 
have a finite group velocity of order ~ A, and are prop- 
agating along the dislocation line, but in the opposite 
directions on each of the dislocation lines. In this sense 
we call each of them a ID chiral mode. Note that here the 
meaning of "chiral" is different from when we used the 
same word to describe the chiral property of the Fermi 
arc surface state or its subbands. Indeed, on each of the 



dislocation line the circular orbital motion of an electron 
in the Fermi arc state (around a hypothetical cylinder 
of radius ro ~ 1) is in the same anti-clockwise direction, 
but following a spiral which evolves in the opposite direc- 
tion (+z or —z) reflecting the opposite direction of the 
Burgers vector. 

VI. CONCLUSIONS 

We have studied electronic states of a 3D Weyl 
scmimetal, which serve, in this regard, as the 3D counter 
part of graphene. Naturally the corresponding Fermi arc 
surface states could be regarded as a 2D version of ID 
edge modes with a flat dispersion, which are known to 
exist in the zigzag edge nano-ribbon. 

This paper, however, points out a crucial difference 
between the two systems. The Fermi arc surface states 
exhibit a specific type of (chiral) spin-to-surface locking. 
This manifests as spin Berry phase when one consid- 
ers a curved surface, e.g., a cylindrical surface. In ID 
edge states of the zigzag nano-ribbon the edge pseudo- 
spin state is determined rather by the structure of the 
edge I22ESI xhe spin Berry phase has been regarded as a 
hallmark property of the helical surface states of a topo- 
logical insulator. In this paper we have demonstrated 
that in the case of 3D Weyl semimetal, the existence of 
peculiar spin Berry phase in the Fermi arc state leads to 
a number of unusual finite size effects: 

1. The nanowire spectrum shows a feature of multiple 
subbands, which is gapped at E = (at the level 
of Weyl points). 

2. In the case of triply periodic surfaceless system, a 
protected gapless chiral mode appears along a dis- 
location line and dominates the low-energy trans- 
port. 

As a general remark, we have emphasized that whenever 
a system bears a surface state involving a spin Berry 
phase, a scries of plaquettes pierced by a 7r-flux tube 
always hosts a zero-energy bound state. The second 
statement above is a specific version of this general phe- 
nomenon in the case of a 3D Weyl semimetal bearing 
Fermi arc states. 
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